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Abstract 

A correlator of the vector current of a heavy quark is computed analyti- 
cally near threshold in the next-to-next-to-leading order in perturbative and 
relativistic expansion that includes a 2 ., a s v and v 2 corrections in the coupling 
constant and velocity of the heavy quark to the nonrelativistic Coulomb ap- 
proximation. Based on this result, the numerical values of the 6-quark pole 
mass and the strong coupling constant are determined from the analysis of sum 
rules for the T system. The next-to-next-to-leading corrections are found to be 
of order of next-to- leading ones. 
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Insufficiency of the ordinary PT for description the near threshold behavior of vacuum 
polarization function was noted long ago in the context of Coulombic resummation in 
nonrelativistic QED [|. Recently a considerable progress has been made in study- 
ing the near threshold production of heavy quark- ant iquark pair within perturbation 
theory of QCD with resummation of threshold singularities. Both perturbative and 
relativistic corrections have been taken into account in the next-to-next-to-leading 
order in the coupling constant and velocity of the heavy quark to the leading non- 
relativistic approximation based on Coulomb potential 0, [|, This theoretical 
development provides more accurate description of the heavy quark vacuum polariza- 
tion function in the threshold region necessary for such applications as the top quark 
production |J and the precise quantitative investigation of the T system J7], |]. In 
the latter case higher order corrections to leading Coulomb behavior in the thresh- 
old region are essential both numerically for extracting the 6-quark mass and the 
strong coupling constant || and qualitatively for justifying the perturbative expansion 
around Coulomb solution. The analytical calculation of the next-to-next-to-leading 
order corrections has not been completed yet though some results are available?]. 

In this paper we present the complete analytical expression for a correlator of the 
vector current of heavy quarks near threshold in the next-to-next-to-leading order 
resumming all 0[(a s /v) n x (a 2 , a s v, v 2 )} terms, with v being the heavy quark velocity. 
The correlator is further used for determination of the bottom quark pole mass m?, 
and the strong coupling constant a s from sum rules for the T system. 

We study the near threshold behavior of the polarization function IT(s) of the 
6-quark vector current j M = 67^6 



within the nonrelativistic expansion || which in the next-to-next-to-leading order 
1 Semi-analytical analysis of the complete next-to-next-to-leading order corrections to the heavy 
quark polarization function near the two-particle threshold has been done in the context of the 




photon mediated t quark pair production (J, ||j 
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reads 



N r ( „ , , „ , , s 4 k 2 



n ( s ) = 7^2 Ca(«-)G(0, 0, k) + ~—G c (0, 0, k) (1) 



with A; = ym 2 — s/A being a natural energy variable near threshold. First term in 
brackets gives the representation for the correlator within NRQCD with Ch{ot s ) being 
a perturbative coefficient matching correlators of relativistic and nonrelativistic vector 
currents. The coefficient Ch{a s ) is computable in full QCD and by now is known to 
the second order in a s expansion 

C h (a s ) = l-C 1 h C F ^ + ClC F (^) 2 

7T V 7T / 

with ^ = 4 and 
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in MS renormalization scheme j|, |5|, [TTJ. Here the group invariants for QCD are 
Ca = 3, C F = 4/3, Tp = 1/2, and 7^ = 0.577216 ... is the Euler constant, Q(z) is the 
Riemann (^-function, rif is the number of light flavors, and /?o = llC^/3 — 4Tpn//3. 
The quantity G(x, y, fc) is the nonrelativistic Green function (GF) of the following 
Schrodinger equation 

/ A x A 2 „ , , a,,, , s /a x 2 
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+V NA (x) + V BF (x, s) + — G(x, y, fc) = 5(x - y) (3) 
where Vc(^) = —Cpa s /x is the Coulomb potential which is supposed to dominate 
the whole QCD interaction in the energy region of interest, x = |x|, Vna( x ) — 
—C aC foi 2 s / {mbX 2 ) is the non-Abelian potential of quark-antiquark interaction ||12|| , 
Vbp(x, s) is the standard Breit-Fermi potential (up to the color factor C F ) containing 
the quark spin operator s, e.g. [I3~|. The terms V% (i — 1, 2) represent first and second 



order perturbative QCD corrections to the Coulomb potential [ 13 , ITS 

V 1 (x) = l^(x)(C 1 + C 1 1 ln(x/i)), 



V 2 {x) = V c {x){Cl + C\ ln(x/i) + C\ ln 2 ^)), (4) 

where 

Cl = a x + 2f3 0lE , C\ = 2/3o, 

Cl ={j + 4 ^) $ + 2 ^ + 2 ^ihE + (h, 
C\ = 2(/3x + 2/3oax) + C 2 = 4/3 2 , 

Ql = ~9~ A ~ ~9 Fnf > 
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The second term in eq. ([!]) is generated by the operator of dimension four in the 
nonrelativistic expansion of the vector current (see, for example, |fL6H ). It contains 
the GF of the pure Coulomb Schrodinger equation W7\ at the origin 



Gc ,0, 0, t) . + „ g) + 7£ + (x - ^5)) ,5) 

where = r'(a;)/r(a;) and T(x) is the Euler T-function. 

The solution to eq. (|3|) can be found within the standard nonrelativistic pertur- 
bation theory around the Coulomb GF Gc(x, y, k). The leading order corrections to 
the Coulomb GF at the origin due to A 2 , Vna and Vbf terms are known analytically 
H HQ. After including these corrections the approximate GF of eq. (H) at the origin 
takes the form Q 



2 The term Vat a can be fully accounted for the Coulomb GF because the corresponding differential 
equation is exactly solvable in standard special functions. Numerically this is not important for 
applications though. 



i 2 b V 2C F , 

where \I/ 2 (^) = ^i( x )- Note that in ref. || the shift of the spectrum of intermediate 
nonrelativistic Coulomb bound states was treated exactly i.e. without expanding of 
the energy denominators. This accounts for a part of the higher order corrections. 
We, however, consistently work in the next-to-next-to-leading order and keep only the 
second order terms in eq. (^j). Since this part of the corrections is relatively small the 
difference between these two approaches is really negligible for the numerical analysis 
of the sum rules. 

The correction AGi to eq. (Q) due to the first iteration of V\ term of the QCD 
potential has been found in ref. where the consistent analysis of sum rules for bb 
system in the next-to-leading order has been performed 

AG 1 (0, 0, k) = (f^ F 2 (m)(m + 1) (c] + (L k + * x (m + 2))C} 

oo m—1 i i oo 

-2 X) E F(m)F(n)—— C\ + 2 £ (C l Q + (L fc - 2 7i; - ^(m + 1))C? 

m=l n =0 m 72 m=0 

+^^0 + (-IeL, + ^L 2 ) C 1 ) (7) 

where L k = In f^M and 

C F a s m b ( C F ct s m h \~ x 
r [m) = — — m + 1 — 



(m + l)2fcV 2A; 
The correction AG® to eq. @ due to V2 part of the potential is also known || 

AGf (0, 0, fc) = (^) 2 (f>V) ((- + 1) (Cg + ^ + L 2 fe Cf) 

+ (m + l)*!(m + 2) (C 2 + 2L fc C 2 ) + J(m)C 2 
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where 
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J(m, n) = 2 *i(m - n) h 27 E 

m — n \ n + 1 
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+2 (^i(m - n + 1) - * x (m + 1)), 

m — n 

K(m) = 2(^i(m + 1) + lE f + ^ 2 (m + 1) - # 2 (m + 1) + 2 7 |, 

In this paper we complete these results by computing the correction AG^ due to 
the second iteration of V\ term which of the proper (next-to-next-to-leading) order 
according to counting of smallness in nonrelativistic QCD with respect to a s and v. 
The result reads 

AG « (0 , ,^ (S) ^ 1 |(£ H 3 (m)(m + 1) 
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(9) 



where 

t_t/ \ / , , C F a s rr, 

H (m) = I m + 1 — 

\ ZK 

We are going to describe the details of this rather cumbersome calculation elsewhere. 
One remark is in order though. Because ultraviolet divergences in eq. depend on k 
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one has to match the calculation of these corrections to the calculation of the Wilson 
coefficient Ch{a B ) (eq. @) 0, |^]. Such matching is not necessary for the calculation 
of AGi and AG% terms because their divergent parts are k independent. 

Thus eqs. fll], 0-|9|) give the complete analytical expressions for the vacuum po- 
larization function of heavy quarks near the two-particle threshold in the next-to- 
next-to-leading ordeif]. Eqs. (||-|9]) look awkward and they can be rendered into more 
readable form by using \1/ functions for expressing some of the sums entering the 
formulae. However for direct numerical analysis of sum rules for bb system this form 
is most suitable with respect to applicability of efficient numerical algorithms of a 
symbolic system. 

Obtained formulae are applied to the analysis of the T system for extraction of the 
6-quark pole mass m& and the coupling constant a s . The sum rules are formulated in 
the literature @, § and we will use the latest version || with correct large n behavior. 
The moments M. n 

R(s)ds 



= (4m 

s=0 



2\n 



s n+1 



of the spectral density R(s) = 127rlmil(s + ie) are compared with experimental ones 

M exp__ rR b (s)ds 



Ql Jo 8"+! 

under the assumption of quark-hadron duality. The experimental moments Ai^" p 
are generated by the function Rb{s) which is the normalizad cross section Rb(s) = 
a(e + e~ — > hadrons b 5)/cr(e + e~ — > Here Q b = —1/3 is the 6-quark electric 

charge. 

Numerical values are obtained basically by saturating the experimental moments 
with the contribution of the first six T resonances (see for details). Their leptonic 



widths Tfc and masses (fc = 1 ... 6) are known with good accuracy [18 



M e XP _ Ml! ( 9?r f -J*- + rds RbisY 



Ql W Q ED{m h )ti M k J* s 



3 In refs. ||, |) the corrections to the Coulomb GF due to Vi (i = 1, 2) terms of the potential were 
treated numerically for complex values of energy far from the real axis. 
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The rest of the spectrum beyond the resonance region for energies larger than so ~ 
(11.2 GeV) 2 (continuum contribution) lies far from threshold and is safely approx- 
imated by the ordinary PT expression for the theoretical spectral density, so there 
Rb{s) ~ R(s). The influence of the continuum on high moments is almost negligible 
numerically and in any case under strict control^. Electromagnetic coupling constant 



is renormalized to the energy of order of mj with the result ctQED( m b) — l-07a 2 [18 

We work with moments for 10 < n < 20 that simultaneously guarantees the small- 
ness of both the continuum contribution and the nonperturbative power corrections 
due to the gluonic condensate ||. The first one is not well known experimentally 
and has to be suppressed to make results independent of sq. The second one should 
be small because the value of gluonic condensate (and higher order condensates) is 
not known well numerically. The normalization point /i = m& is used throughout the 
computation^]. At the scale \i = nib both the hard and soft gluon corrections are of 
the same order of magnitude. For a lower scale the hard corrections become large 
while for higher scale the same is true for the soft corrections ||. We found that at 
H ~ nib the \i dependence of the results is minimal which is a solid indication that at 
this point the higher order corrections are also small. 
The result of the fit is 

a s (m b ) = 0.22 ±0.02, or a s (M z ) = 0.118 ± 0.006. 

The sum rules are much more sensitive to the 6-quark mass than to the strong coupling 
constant so it is instructive to fix a s [Mz) = 0.118 to the "world average" value JT8| 



The expressions for the first few moments of the spectral density are now available in ordinary 

perturbation theory with a 2 s accuracy fl9|| , however, they cannot be used in theoretical formulas 

for sum rules directly because the spectrum is well known experimentally only for energies close to 

threshold due to existence of sharp resonances while the contribution of the continuum to these low 

moments is large in comparison with the resonance contribution. 

We work strictly in the next-to-next-to-leading order approximation and, therefore, use the same 

normalization point for soft and hard corrections in contrast to [Q, || where different normalization 

points were chosen for these two parts. 
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and then to extract m b {n). In this way we obtain the following estimate for the mean 
value over the considered range of n 

m b = 4.78 GeV. 

This value is in a good agreement with the results of the first order analysis || where 
at a s {Mz) = 0.118 we obtained 

m b = 4.75 GeV. 

Note that the optimization procedure |?(| was used to improve convergence of per- 
turbation theory in the previous analysis ||. As we see this procedure turns out 
to be a powerful tool to estimate the higher order contributions. For comparison, 
the leading order result is m b = 4.70 GeV and in the next-to-leading approximation 
without optimization one gets m& = 4.72 GeV. 

Main uncertainties of numerical values for considered parameters stem from the 
same sources that were identified in ref. 0. The error coming from n distribution 
for the mass at fixed value of the coupling constant is about ±0.5% for 10 < n < 20. 
The \x dependence for /i = m b ± 1 GeV (where this dependence is minimal) introduces 
another ±0.5% of uncertainty. Thus our final estimate of the bottom quark pole mass 
is 

m b = 4.78 ± 0.04 GeV. 

Note that the uncertainty originated from the n and \i dependence is not reduced 
in comparison with the next-to-leading order. This means that the contribution of 
the higher order corrections which has to cancel n and \x dependence of the results 
is still important. Let us emphasize that the convergence of the perturbation the- 
ory for the vacuum polarization function of heavy quark near threshold is not fast. 
We have found the next-to-next-to-leading order corrections to be of order of the 
next-to-leading ones. Furthermore, in the case of 6-quark the corrections due to 
the perturbative modification of the Coulomb instantaneous potential (i.e. related 
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to AGi and AG2 terms) dominate the total correction in the next-to-leading and 
next-to-next-to-leading orders. Inclusion of these corrections is quite important for 
consistent analysis of sum rules for the T system. 

To conclude we have constructed an expression for the vacuum polarization func- 
tion of the vector current of a heavy quark near threshold. It is completely analytic 
in the next-to-next-to-leading order in perturbative and relativistic expansion up to 
a 2 , a s v and v 2 corrections. The polarization function was used for determination of 
the 6-quark pole mass and the coupling constant from sum rules for the T system 
that are saturated by contributions near threshold. In fact, there is no much hope for 
improving our results: next order approximation seems to be too complicated for an- 
alytical treatment within the regular perturbation theory for NRQCD. The analysis 
showed a remarkable stability with respect to the next-to-leading one supplied with 
an optimization procedure in a variational spirit. Having in mind the considerable 
technical difficulty of computing next approximation and recognizing the necessity of 
improving the theoretical predictions in view of new high quality experimental data 
we think that the next step in the near future will be connected with optimization of 
the present approximation. 

Acknowledgements 

We thank J.H.Kuhn for support, encouragement, and discussions. A.A.Penin great- 
fully acknowledges discussions with K.Melnikov. This work is partially supported by 
Volkswagen Foundation under contract No. 1/73611. A.A.Pivovarov is supported in 
part by the Russian Fund for Basic Research under contracts Nos. 96-01-01860 and 
97-02-17065. The work of A.A.Penin is supported in part by the Russian Fund for 
Basic Research under contract 97-02-17065. 



References 

[1] M.A.Braun, ZhETP Lett. 27(1968)652. 



10 



[2] R.Barbieri, P.Christillin and E.Remiddi, Phys.Rev. A8(1973)2266. 

[3] J.H.Kiihn, A.A.Penin and A.A.Pivovarov, Preprint TTP-98-01, 
|hep-ph/9801356 . 



[4] A.H.Hoang and T.Teubner, Preprint UCSD/PTH 98-01, [hep-ph/ 980 1397 



[5] K.Melnikov and A.Yelkhovsky, Preprint TTP-98-10, frep-ph/9802379 . 



[6] V.S.Fadin and VA.Khoze, Pis'ma Zh.Eksp.Teor.Fiz 46(1987)417; Yad.Fiz. 
48(1988)487; 

W.Kwong, Phys.Rev. D43(1991)1488; 

M.J.Strassler and M.E.Peskin, Phys.Rev. D43(1991)1500; 

M.Jezabek, J.H.Kiihn and T.Teubner, Z.Phys. C56(1992)653; 

Y.Sumino, K.Fujii, K.Hagivara, H. Murayama and C.-K.Ng, Phys.Rev. 

D47(1993)56; 

K.Fujii, T.Matsui and Y.Sumino, Phys.Rev. D50(1994)4341. 

[7] V.A.Novikov et al, Phys.Rev. Lett. 38(1977)626; 
V.A.Novikov et al, Phys.Rep. C41(1978)l; 
M.B.Voloshin, Yad.Fiz. 36(1982)247; 

M.B.Voloshin and Yu.M.Zaitsev, Usp.Fiz.Nauk 152(1987)361. 

[8] M.Voloshin, Int.J.Mod.Phys. A10(1995)2865. 

[9] W.E.Caswell and G.E.Lepage, Phys.Lett. B167(1986)437. 

[10] G.Kallen and A.Sarby, K.Dan. Vidensk.Selsk.Mat.-Fis.Medd. 29(1955), N17, 1. 

[11] A.H.Hoang, Phys.Rev. D56(1997)7276; 

A.H.Hoang, J.H.Kiihn and T.Teubner, Nucl.Phys. B452(1995)173; 
A.Czarnecky and K.Melnikov, Preprint TTP-97-54, |hep-ph/971222^ ; 



M.Beneke, A.Signer and V.A.Smirnov, Preprint CERN-TH-97-353, 
hep-ph/9712302 . 



11 



[12] S.N.Gupta and S.F.Radford, Phys.Rev. D24(1981)2309; Phys.Rev. 
D25(1982)3430 (Erratum); 

S.N.Gupta, S.F.Radford and W.W.Repko, Phys.Rev. D26(1982)3305. 

[13] L.D. Landau and E.M.Lifshitz, Relativistic Quantum Theory, Part 1 (Pergamon, 
Oxford, 1974). 

[14] W.Fisher, Nucl.Phys B129(1977)157; 
A.Billoire, Phys.Lett. B92(1980)343. 



[15] M.Peter, Phys.Rev.Lett. 78(1997)602; Preprint TTP-97-03, [hep-ph/9702245 

[16] G.T.Bodwin, E.Braaten and G.P.Lepage, Phys.Rev. D51(1995)1125. 

[17] J.Schwinger, J.Math.Phys. 5(1964)1606. 

[18] Particle Data Groop, Phys.Rev. D54(1996)l. 

[19] K.G.Chetyrkin, J.H.Kuhn and M.Steinhauser, Phys.Lett. B371(1996)93; 
Nucl.Phys. B482(1996)213. 

[20] A.A.Penin and A.A.Pivovarov, Phys.Lett. B367(1996)342. 



12 



